Dynamic light diffusion, Anderson localization and lasing in disordered inverted opals: 

3D ab-initio Maxwell-Bloch computation 
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We report on 3D time-domain parallel simulations of Anderson localization of light in inverted dis- 
ordered opals displaying a complete photonic band-gap. We investigate dynamic diffusion processes 
induced by femtosecond laser excitations, calculate the diffusion constant and the decay-time dis- 
tribution versus the strength of the disorder. We report evidence of the transition from delocalized 
Bloch oscillations to strongly localized resonances in self-starting laser processes. 



The first investigation of strong localization of light 
in periodically structured 3D materials dates back to 
more than twenty years ago [l|. Since then the field 
of photonic crystals (PhC) has steadily grown, and 
evolved into various directions including all-optical sig- 
nal processing, quantum information and bio-sensing 
@, H H S, @, 0, @, 0. Notwithstanding these develop- 
ments, the issue of 3D light localization in the presence of 
disorder and backbone periodicity is still largely debated. 
As randomness is unavoidable in any practical realization 
of PhC, this topic has a fundamental importance for the 
physics and the applications of these devices. 

The starting point is a medium displaying a complete 
photonic band-gap (PBG), i.e. a spectral region where 
the density of electromagnetic states vanishes. If disorder 
is introduced, the PBG is progressively filled by tails of 
resonances at its edges 1(| 11. 12|; this circumstance nur- 
tures the observation of 3D Anderson (or strong) local- 
ization of light, which has demonstrated to be elusive in 
absence of periodicity 0, Q EH, H Ezl 01 • Indeed the 
PBG is filled by states that are expected to be exponen- 
tially localized inside the PhC; following [H, these should 
induce a diffusive regime of light propagation 19|, 13, 2l[ 
and a "critical slowing down of photons", which results 
in the development of an exponential temporal tail for a 
pulsed beam propagating in the PhC. This effect, known 
as dynamic light diffusion, is currently largel y in vesti- 
gated in disordered systems [13, E3, E3, 0, 0. 

For what concerns 3D PhC, however, there are several 
open problems that are extremely relevant for the final 
assessment of the observation of Anderson localization of 
light that is still missing. Open fundamental questions in- 
clude the possibility of exciting these resonances by short 
pulses and, correspondingly, the expected values of mea- 
surable quantities such the dynamical diffusion constant 
and the decay-time distribution of the transmission (i.e. 
the way the latter deviates from a simple exponential). 

A key issue is the role of disorder-mediated resonances 
in self-starting lasing (see, e.g, [27]): if a small per- 
centage of randomness can induce their existence, any 
practical light-emitting device based on 3D PhC (as in 
[281 . [IK) is expected to be strongly affected by Ander- 



son localizations. Therefore one can argue about their 
role in the stimulated emission processes if the PhC is 
filled by a light- amplifying material. Are 3D Anderson 
localizations able to sustain a laser action? Or delocal- 
ized "Bloch" modes always prevail? From a theoretical 
perspective, the current state of the art is strongly ham- 
pered by the huge amount of the needed computational 
resources when dealing with 3D problems. Indeed sim- 
pler low-dimensional models (in ID and 2D) are quite 
unsatisfactory for many important reasons, including the 
lack of a complete PBG and the critical nature of Ander- 
son localizations in 3D [301 ]. 

In this manuscript, we apply a massively parallel ap- 
proach to give answers to the mentioned open problems. 
For a realistic finite structure, we show the way the PBG 
is closed when increasing the strength of the disorder, we 
provide evidence of the excitation of localized modes by 
using short pulses, and we determine the relation between 
the diffusion constant D and the strength of the disorder. 
We find that, in the strongly localized regimes (when D 
vanishes) , the decay-time distribution of the transmission 
splits into a fast and a slow component, showing a clear 
dynamical signature of the existence of a second regime of 
light localization in disordered PhC. Finally, in order to 
ascertain if Anderson resonances are actually involved in 
lasing processes, we solve from first principles the most 
general semi-classical model of light-matter interaction 
to simulate a self-starting laser oscillation. Specifically, 
we show the transition from emission mediated by delo- 
calized periodical waves (Bloch modes) to that driven by 
Anderson localization. 

Sample and numerical approach — Our approach 
is based on a parallel Finite-Difference Time-Domain 
(FDTD) code [3jJ with typical runs employing from tens 
to hundreds of processors on Linux clusters and IBM-SP5 
systems. With reference to lasing processes, we solve the 
Maxwell-Bloch equations for a four-level system within 
a first-principle formulation (32[. At variance with pre- 
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viously reported ID and 2D models | 
levels are taken as degenerated and coupled to the three 
components of the electric field, thus accounting for a 
fully vectorial representation of the light interaction with 



2 



a resonant system. The density matrix equations for a 
four level systems are re-written in terms of real- valued 
fields by applying an homomorphism with the su(4) al- 
gebra [32[ and solved in the middle of the Yee grid. The 
parameters are chosen in order to mimic a typical dye, 
i.e. Rhodamine B, with resonance wavelength 484THz 
(Ao = 620nm) and refractive index n = 1.3. The charac- 
teristic time-scales for the density matrix equations are 
set to T\ — 1 ns and T2 = 10 fs, providing a small 
signal gain bandwidth of the order of I/T2 = lOOTHz. 
The noise due to spontaneous emission is included as a 
stochastic source term in the electric field evolution [35[ . 
The integration domain is a box with size L^ ox = 6/im 
and "UPML" boundary conditions [HI. Most of the re- 
ported simulations adopt grid sizes of 200 3 , with 400 3 
high-resolution runs for checking. The considered finite- 
size PhC (edge L = 4/im) is an inverted opal (Fig Q] inset) 
with dielectric constant is e r = 12, displaying a complete 
PBG around a/A = 0.8 where a is the cubic lattice con- 
stant. In absence of disorder, the radius of each sphere 
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in the PhC is r = a/2 a/2 [36 
suits are reported with reference to the case a = 480nm, 
providing a band gap around 500THz (A = 600nm). Dis- 
order is finally introduced by a random-number genera- 
tor providing uniform deviates £ in [—1/2, 1/2], and such 
that the radius of each air sphere in the inverted opal is 
given by r = 7*0 (1 + 7£), where 7 measures the strength 
of the disorder (7 = is the ordered case). 

Wide-band excitation and gap filling. In the first set of 
runs we consider the passive air-sphere PhC and investi- 
gate the way the forbidden band is closed when increas- 
ing the strength of the disorder. We launch a single cycle 
(linearly polarized TEMqq with waist 2 jam) pulse with 
frequency 484THz (A = 617.25nra) in the [100] crystal- 
lographic direction (denoted as z hereafter), and we re- 
trieve the field in a low symmetry point inside the sample. 
After a total time interval of 2.9ps (10 5 time-steps with 
200 3 grid and dt = 2.9 x 10~ 17 s), we calculate the electro- 
magnetic oscillation spectrum with resolution ~ 0.4 THz. 
Figure [T] shows the field spectrum for various 7, with re- 
sults averaged over 10 realizations of the disorder (the 
total processor time for figure [T] is of the order of 480 
days). For large 7 the PBG is closed by tails of states at 
the band-edges. This demonstrates that of a few percent 
of randomness is sufficient to close the forbidden band. 

Dynamical diffusion and localization. In a second set 
of runs we consider the dynamic light diffusion process of 
a 100/5 Gaussian pulse with spectrum located within the 
PBG in proximity of the the high-frequency edge (FigHJ 
other beam parameters are the same as the single-cycle 
case considered above). We record the transmitted power 
pulse T(t) (transversely integrated Poynting vector at the 
output face of the PhC as in [26[). Figure [2^ shows T(t) 
for increasing 7: we observe the development of a pro- 
nounced temporal tail, which unveils the "critical slowing 
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FIG. 1: (Color online) Spectral content (log. scale, each 
line is vertically shifted by an arbitrary amount) of the elec- 
tromagnetic field inside the 3D inverted opal for increasing 
degree of disorder; the thick line is the ordered case (7 = 0). 
The dot-dashed line is the spectrum (arbitrary linear scale) of 
the pulsed lOOfs excitation for the dynamic diffusion analysis 
(see text). Inset: sketch of the 3D PhC index profile. 
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down of photons" [II] ; this is also evident in the inset [2^ 



FIG. 2: (Color online) (a) Tail of the transmitted pulse for 
various 7 (semi- log. scale, dashed line for 7 = 0), the time 
axis has been shifted such that the peak of T(t) corresponds 
to t = 0; the inset shows the transmitted pulse for 7 = 0.1 
(light region, total time window 2ps) with the input pulse 
superimposed (dark region); (b) Snapshot of E y component 
in the PhC middle (y, z) section for 7 = and t = 6ps; (c) as 
in (b) for 7 = 0.1. 



by comparing the input pulse with the transmitted one 
in the case 7 = 0.1. Figure [2)3, c show the snapshot of 
the E y component of the electric field in the middle as- 
sertion plane (y, z) taken at very large times from the 
peak value of the input pulse, (in this way we select the 
long living modes in the structure). Note that in the ab- 
sence of disorder, (Fig|2}3) a delocalized periodical Bloch 
mode lying in at the band-edge survives; it corresponds 
to the high-frequency spectral tails of the input pulse. 
Time-domain spectral analysis (not reported) furnishes 
its frequency / = 490THz (A = 614nra). In the presence 
of disorder 7 = 0.1 (FigEt), the field snapshot at long 
times reveals the onset of a localized states [1] that turns 
out to correspond to / ^ 484THz (A 9* 620nm). 

Diffusion constant and decay times. Figure [3fci shows 
the dynamic light diffusion constant D versus the degree 
of disorder. We calculate D by fitting the trailing edge 



FIG. 3: (Color online) (a) Diffusion constant D (dots, left- 
scale), quadratic polynomial best-fit (continuous line), and 
time delay To (squares and best-fit dashed line, right scale), 
(b) Decay-time distribution g(r) for various 7, all the lines 
have been vertically scaled to the peak value. 

of T(t) by an exponential tail exp(— t/ro), where To is 
the leading time constant; D is then calculated by D = 
L 2 /7r 2 To. A quadratic fit for D gives 

D^s' 1 ) ^ 14 - 197 - 204 7 2 . (1) 

Equation (pQ) implies that D = when 7 = 0.2. At that 
value for 7 the PBG is closed (FigH]) and a simple expo- 
nential fit for the tail of T(t) does not work. We study 
this regime by the distribution g(r) of the decay-times 
(FiglSb). We find g(r) by minimizing the difference be- 
tween the trailing edge of the transmission T(t) and a 
signal given by the sum of exponential functions with 
decay-constant r and weight g(r). While increasing 7, 
g(r) spreads around the leading diffusion time To- There- 
fore, T(t) progressively deviates from a simple exponen- 
tial unveiling the excitation of localized states. Notably 
enough, as 7 > 0.15 the time constant distribution is 
splits into a fast and a slow component (FigJSb). This 
behavior is due to the fact that when the PBG is closed, 
some of the input energy is trapped in the localized states 
[thus providing the long time tails in T(t)], while other 
travels by tunneling through low-Q factor modes, which 
do not exist at small 7. Previous investigations on low- 
index contrast PhC [42[ envisaged the existence of a sec- 
ond localization regime in a complete PBG in the pres- 
ence of large disorder. Our analysis shows that there is 
a transition in the decay-time distribution in correspon- 
dence of the zero of the diffusion D, when the PBG is 
practically filled by tails of localized states. This is an 
experimentally testable signature of a specific regime of 
light localization, which appears as a characteristic fea- 
ture of PhC. 

3+ ID simulation of lasing in disordered and ordered 
PhC To model laser emission in the PhC, we take the air- 
spheres as filled by an active resonant medium (modeled 
by the Maxwell-Bloch equations detailed above). We first 
analyze the way the presence of the infiltrated medium 
affects the width of the forbidden-gap. Figured! compares 
the responses of an air-filled (dashed-line) and dye-filled 
(thick-line) opal to a single-cycle pulse. The forbidden- 
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FIG. 4: (Color online) (a) Wide-band single-cycle excitation 
spectrum (log. scale) for 7 = 0: filled (continuous line) and 
the empty opal (dashed line). The spectra have been verti- 
cally translated of an arbitrary amount for the sake of compar- 
ison, (b) Self-starting laser spectrum in the ordered (dashed 
thin line) and in the disordered case (thick line) (7 = 0.1). 

band is progressively closed as the refractive index of 
the infiltrated medium increases; indeed we found that a 
"small" gap survives around / = 480THz (A = 625nm). 
We then consider the case the active medium is externally 
pumped. At low inversion 7V a , the laser is below thresh- 
old and only fluorescence is obtained (not reported). At 
high pumping (N a = 5 x 10 25 m -3 ), in the ordered case 
(7 = 0), we find that the PhC cavity starts oscillating 
on a spectrum mainly concentrated at the band-edges of 
the linear PBG of the infiltrated opal (see Fig. SJd). This 
spectrum corresponds to the excitation of Bloch-modes, 
whose spatio-temporal profile is delocalized. This is con- 
firmed by FigEK, showing a 3D iso-surface plot of the 
electromagnetic energy density stored inside the struc- 
ture at long times (t = 3 ps) since the turning on of 
the pump. Conversely, in the disordered case (7 = 0.1), 
we find that the emission spectra is mostly concentrated 
within the linear band gap of the filled case (Fig. [4]). The 
corresponding mode profile is clearly localized, as shown 
in the panel (b) in Fig|5l 

Conclusions — We reported on an extensive investi- 
gation of the interplay between localization and disorder 
in an inverted opal exhibiting a complete photonic band 
gap. The dynamic diffusion regime in response to an 
in-gap lOOfs light pulse excitation shows that the dif- 
fusion constant decreases when increasing the strength 
of the disordered and localized modes are put into os- 
cillation. The time-dependent transmission develops a 
non-exponential tail, evidencing of a second regime char- 
acterized by the splitting of the decay-time distribu- 
tion at large disorder. The fully vectorial 3D Maxwell- 
Bloch simulation of self-starting lasing processes finally 
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FIG. 5: (Color online) 3D iso-surface for the electromagnetic 
energy density at high pumping rate: (a) ordered 7 = and 
(b) disordered case 7 = 0.1. Opaque surface corresponds 
to the energy density level of 3fJ//im 3 ; transparent surface 
corresponds to 0.2fJ//im 3 ; the lateral plots are projection of 
the middle sections of the energy density. 



demonstrates that, in the presence of disorder, delocal- 
ized Bloch-modes are replaced by the Anderson localiza- 
tions in sustaining stimulated emission processes. 
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